On the multifractal statistics of the local order parameter at random critical points 

application to wetting transitions with disorder 



Cecile Monthus and Thomas Garel 
Servtce de Physique Theorique, CEA/DSM/SPhT 
Unite de recherche associee au CNRS 
91191 Gif-sur-Yvette cedex, France 

Disordered systems present multifractal properties at criticality. In particular, as discovered by 
Ludwig (A.W.W. Ludwig, Nucl. Phys. B 330, 639 (1990)) on the case of diluted two-dimensional 
^— ■») ' Potts model, the moments p'^ir) of the local order parameter p{r) scale with a set x[q) of non-trivial 

, exponents x{q) 7^ qx{l). In this paper, we revisit these ideas to incorporate more recent findings: 

■ (i) whenever a multifractal measure w(r) normalized over space '^^'wir) — 1 occurs in a random 

system, it is crucial to distinguish between the typical values and the disorder averaged values of 
the generalized moments Yq — X]r^''('')' since they may scale with different generalized dimensions 
D(q) and D(q) (ii) as discovered by Wiseman and Domany (S. Wiseman and E. Domany, Phys Rev 
E 52, 3469 (1995)), the presence of an infinite correlation length induces a lack of self- averaging 
at critical points for thermodynamic observables, in particular for the order parameter. After this 
general discussion valid for any random critical point, we apply these ideas to random polymer 
models that can be studied numerically for large sizes and good statistics over the samples. We 
I study the bidimensional wetting or the Poland-Scheraga DNA model with loop exponent c — 1.5 

^ , (marginal disorder) and c = 1.75 (relevant disorder). Finally, we argue that the presence of finite 

Griffiths ordered clusters at criticality determines the asymptotic value x{q — > 00) — d and the 
minimal value Ofmin = D{q 00) = d — x{l) of the typical multifractal spectrum /(«). 
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I. INTRODUCTION 



Among the various areas where multifractality occurs (see for instance [H, [2, H, 0, [E S and references therein), 
P5 the case of critical points in the presence of frozen disorder is of particular interest. The idea that multifractality 
O occurs at criticality has been first established for quantum Anderson localization transitions [1, Q and has been the 
1^1 ' subject of very detailed studies [Tol. [Til. [T^ . For the directed polymer in a random medium in dimension 1-1-3, where a 
[ disorder-induced localization/delocalization occurs, the multifractal properties studied recently [3] are very similar to 
T— I ■ the case of Anderson transitions. In the field of spin models, the most studied case seems to be the two-dimensional 
^ ' diluted g-state Potts model, where multifractality was discovered by Ludwig [13] via conformal field theory using 
. perturbation theory in the parameter (o — 2) governing the disorder relevance. This work has motivated numerical 
' studies for various values of q [l^, [l^, [13, EBl ■ The idea of multifractality has been also proposed in other disordered 
. models like spin-glasses and random field spin systems fl^, [1^ |2l[ , and has been studied numerically for spin glasses 
• ' on diamond hierarchical lattice ^22] . Finally, for disordered quantum spin-chains, it turns out that the statistics of 
1^ I critical correlation functions is described by "multiscaling" , which is even stronger than multifractality l23l . This is 
' because these disordered quantum spin-chains are actually governed by "Infinite disorder fixed points" |24| . 
I So the presence of multifractality at criticality seems generic in disordered systems. However, following p^ . 
• • . most studies on classical disordered models have focused on the statistics of two-point correlation functions, whereas 
multifractality already occurs at the level of one-point functions like the order parameter or the energy density [T3 | . In 
particular, the moments of the local order parameter scale with a set x(n) of non-trivial exponents x(n) ^ nx(X). In 
^ ■ this paper, we revisit these ideas in the light of more recent findings, concerning the possible differences between the 
5^ 1 exponents for typical and averaged values, and the lack of self-averaging of thermodynamic observables at criticality 
[25I [2^ W\ ■ We then study the multifractal statistics of the order parameter in random polymer models that can be 
studied numerically for large sizes and good statistics over the samples. 

The paper is organized as follows. In Section |TT1 we summarize the outcome of previous works concerning the 
multifractal statistics and the lack of self-averaging at random critical points. In the remainder of the paper, we apply 
these ideas to wetting transitions with disorder. These random polymer models are presented in Section IIIII Our 
numerical results on the statistics of the local order parameter are given respectively in Sections IIVI for loop exponent 
c = 1.75 (relevant disorder) and in Section |V] for loop exponent c = 1.5 (marginal disorder). In Section IVll we discuss 
the influence of boundary conditions on the multifractal spectrum. Section [VIII contains our conclusions. 
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II. MULTIFRACTAL STATISTICS OF THE LOCAL ORDER PARAMETER AT CRITICALITY 
A. Disorder averaged moments of local order parameter 

Let p{r; i, L) be the local order parameter at site r, in a finite disordered sample i of volume L"* in dimension d. 
This local order parameter is usually defined in terms of a thermal average, for instance p{r;i,L) —< a{r;i,L) > 
in disordered ferromagnets, p{r]i,L) —< a{r\i,L) >^ in spin-glasses. For the random polymer models described in 
Section imi p(r; i, L) corresponds to the contact density of monomer r (Eq. H^ . Let L) denote the spatial average 
over all points r of the sample 

In a pure system, the exponent Xpure that governs the decay of the spatial average Ppure{L) also describes the decay 
of the local parameter Pp„re(^; L) for any point r in the bulk 

Ppure{L) ~ ~ Ppure{r G bulk; L) (2) 

In a disordered sample however, the spatial heterogeneity of the disorder induces a spatial heterogeneity in the local 
order parameter p{r] i, L) at criticality. In particular, there exists a family of non-trivial exponents x{q) ^ qx{\) [T^ 
for the disorder averaged powers of the local order parameter 



[p{r;i,L)]i oc — ^ (3) 

In this formula, it is convenient to consider q as a continuous real parameter to probe also non-integer moments and 
even negative moments. 



B. Introduction of a normalized multifractal measure 

Since the multifractal formalism is usually defined for a normalized probability measure llj, it is convenient to 
construct a probability measure from the non- normalized observables one is interested in (ToUlll. [28l [29j . Here for 
the local order parameter, one defines in each sample {i) the following spatial weights 

normalized to 

^u;(r;z,L) = l (5) 

So w(r; i, L) represents the contribution of the site r to the order parameter of the sample i of size L"^. The statistics 
of these weights can be studied via the following generalized moments 

Y,{^,L)^ J2 [w{r■,^,L)Y (6) 

Deep in the ordered phase, where the order parameter p{i, L) is finite as L — > oo, the weights are expected to be of 
the same order (see Eq. 0]). The decay of the generalized moments then follow the simple scaling 

Yq{i, L)\orderedphase ^ ]^(q—l)d ^'^^ 



At criticality however, the generalized moments Yq will display multifractality, with a priori different exponents for 
typical and averaged values. 
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C. Typical generalized dimensions D{q) 

At criticality, the decay of typical values define a series of generalized exponents T{q) — {q — l)D{q) 

Y*yP(L) = e'"^''(''-^) ^ = I (8) 

The exponents D{q) represent generalized dimensions [H : -D(O) represent the dimension of the support of the measure, 
here it is simply given by the space dimension 

D(0) = d (9) 
D{1) is usually called the information dimension [1] , since it describes the behavior of the 'information' entropy 

s{i,L) = -J2 'w{r;i,L)\nw{r;i,L) = --dqYg{i, L)\q^i ~ i:>(l)lnL (10) 

Finally D{2) is called the correlation dimension [l| and describes the decay of 



r2*^P(L) EEe''^^^^^*'^) (11) 

D. Typical singularity spectrum f{a) 

In the multifractal formalism, the singularity spectrum f{a) is given by the Legendre transform of T{q) ^ via the 
standard formula 

q =/'(«) (12) 

riq) =aq-f{a) (13) 

The physical meaning of f{a) is that the number Ml{oi) of points r where the weight w(r; i, L) scales as L~" typically 
behaves as 



A/'*^^'(a) = e'"^^(") oc (14) 
So the Legendre transform of Eq. (fT5)) corresponds to the saddle-point calculus in a of the following expression 

Y^yP{L) ^ ( da L-'" (15) 



The general properties of the singularity spectrum /(a) are as follows [l| : it is positive f{a) > on an interval 
[oi-mim Oimax] whcrc amin ~ D{q — +oo) is the minimal singularity exponent and amax = D{q — — oo) is the maximal 

singularity exponent. It is concave /"(a) < 0. It has a single maximum at some value ao where f{ao) = D{q = 0), 
so here (Eq. [Q]) 

f{aa) = D{q = Q)=d (16) 
The singularity exponent ao is thus the typical value 

"0 = atyp (17) 

However, the singularity that yields the leading contribution to the normalization Yi — w{r) = 1 of the measure 
is the singularity exponent given by the information dimension D{\) of Eq. [TO] 

ai^f{ai)^D{l) (18) 
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E. Generalized dimensions D{q) defined from disorder averaged values 

Following [H, many authors consider that the singularity spectrum has a meaning only for /(a) > (lol . 
However, when multifractality arises in random systems, disorder-averaged values may involve other generalized 
exponents [s^, [Ml, [H, HI] than the typical values (see Eq. [5]). In quantum localization transitions, these exponents 
were denoted by f{q) = {q — l)D{q) in and we will follow these notations 



For these disorder averaged values, the corresponding singularity spectrum /(a) defined by 

JJliaj oc (20) 
may become negative /(a) < [13, [13, HU, [H, [s^l to describe rare events. 

F. Wiseman-Domany lack of self-averaging at criticality 

To make the link between the exponents x{q) defined from the powers of the local order parameter (Eq. [3]) and 
the multifractal exponents of the normalized weights (Eq. 3]) [281], one needs to use the equivalence between spatial 
average and disorder averages for the local order parameter. 

1. Lack of self- averaging of extensive thermodynamic observables 

In disordered systems off-criticality, the densities of extensive thermodynamic observables are self-averaging, because 
the finiteness of the correlation length £,{T) allows to divide a large sample into independent large sub-samples. At 
criticality however, this 'subdivision' argument breaks down because of the diver gen ce of C(Tc) = oo at Tc, and a lack 
of self-averaging has been found at criticality whenever disorder is relevant [25l |26|. \2T\ . More precisely, for a given 
observable X, it is convenient to define its normalized width as 



RxiT,L)^^^^m=mE^ (21) 

(X.(T,i))2 

To be more specific, in ferromagnets, the observable X can be the magnetization M, the susceptibility %, the singular 
parts of the energy or of the specific heat [23]. In terms of the correlation length ^(T), the following behaviour of 
Rx{T,L) is expected : 

(i) off criticality, the correlation length £,{T) is finite. For L ^ £,{T), the system can be then divided into nearly 
independent sub-samples and this leads to 'Strong Self- Averaging' 

Rx (T, £) - ^7 off criticality for L > ^(T) (22) 

(ii) in the critical region, when L ^ £,iT), the system cannot be divided anymore into nearly independent sub- 
samples. In particular at Tc where ^(Tc) = oo, one can have either 'weak self-averaging' for irrelevant disorder 
according to the Harris criterion [4l| , i.e. whenever the pure specific heat exponent apure — 2 — di^pure is negative 

Rx {Tc{oo) , L) ^ L ""p"^' for irrelevant disorder (ttpure < 0) (23) 

or 'No Self- Averaging' 

Rx{Tc{oo), L) ^ Cst for random critical points (24) 

Note that for the marginal case apure = from the point of view of the Harris criterion, the power governing the 
'weak self-averaging' of Eq. E51 vanishes, so the ratio Rx(Tc{oo),L) can either remain finite as in Eq. [M]or vanish 
logarithmically. 
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2. Application to the powers of the local order parameter 

Let us now apply these results to the spatial averages of powers of the local order parameter 

p,(^,L)^-l ^ [p{r;^,L)r (25) 

reL'' 

that generalizes Eq. [l]to arbitrary q. 

In the ordered phase, the 'Strong Self- Averaging' property of Eq. means 

p^i^,L)\T<T.^^^r,iT) + ^ (26) 

where the leading term rq(T) is non-random and coincides with the disorder-averaged value in the thermodynamic 
limit L — > oo 

r,(T) = Hm ([p(r; z, L)]*) (27) 

and where Vq{i) is a random variable depending on the sample (i). 

At criticality, the 'No Self- Averaging' result of Eq. [M] means that the spatial averages defined in Eq. [55] behave 
asymptotically as 



where the exponent x{q) is the exponent governing the decay of the disorder-averaged q-moment p'^(r; i, L) of Eq. [31 
and where Uq{i) is a random variable of order 0(1) depending on the sample (i). 

In the following, we will use these result to understand the relations between the exponents for non-normalized 
observables and for the normalized measure. It will be useful to introduce the rescaled variable 

M,(i,L) =i"(«Vg(«,^) (29) 
that remains a random variable Uq{i) of order 0(1) in the limit L oo. 



G. Relation between the exponents x{q) and r(g) 
In terms of the local order parameter p{r; i, L), the generalized moments Yq(i, L) reads from Eqs[3|and[ni 

Yq{^, L) = ^ (30) 




From Eq [29] concerning the spatial averages of Eq. [25l one obtains 

' ' (L'^-(i)«i(z, L))' " iu^i^, L)Y ^'^^ 

The typical values of the random variables Uq{i, L) in the limit L oo are of order 0(1) and thus the exponents T(q) 
governing the typical values of Eq. [Sj read 

Tiq) = {q- l)D{q) = x{q) - d + q{d - x{l)) (32) 

or equivalently the typical generalized dimensions read 

,,(,)^d„^^im£_^ (33) 

So these relations written in Ref [1^ relate the exponents x{q) of disorder-averaged moments of the local order 
parameter (Eq. [3| ) to the typical exponents T{q) of the normalized measure (Eq. [5]). However, the exponents T{q) 
of Eq. [T9l cannot be simply related to {x{q),T{q)), since the disorder average of Eq. [31] may involve a L-dependent 

rare event contribution of the random variables j^^fiJY^^ ™ particular for large q, since q enters as a power in the 
denominator. 



6 



H. Conclusion for the statistics of the local order parameter 

Let us now come back to our starting point, namely the local order parameter p(r; z, L) at site r. Using Eq. 1291 for 
(7=1, we obtain in terms of the weight of Eq. 0] 

p{r\ i, L) = w{r; i, L) [L"^ p{i, L)] = w{r; i, L)L'^-''^^^ui{i, L) (34) 

where ui(i,L) is a random variable of order 0(1). So the interpretation of the singularity spectrum f{a) for the 
weights given in Eq. [14] can be rephrased as follows : Ml{oi) oc L-^(") represents the the number of points r where the 
weight w{r; i, L) scales as i.e. the number of points r where the local order parameter scales as p(r; i, L) ~ L~y 
with 

y = a-d + x{l) (35) 
In particular, the typical exponent ytyp governing the logarithmic average 

\np{r;i,L) oc -ytyplnL (36) 

L — ^CJO 

is related to the typical value (Xtyp of Eq. [17] by 

ytyp = atyp - d + x{l) (37) 

Similarly, the minimal ymin and maximal ymax exponents are related to amin — D{q =^ +00) and amax — D{<1 = —00). 
In particular, since the minimum value ymin cannot be negative, one has the bound 

ymin = aynm - d + x{l) > Q (38) 



I. Critical region 



Both in quantum localization [lO|, [Uj, [12] and in disordered ferromagnets [IJI , the multifractal statistics exactly at 
Tc is expected to coexist with a single correlation length exponent v outside Tc- More precisely, the powers of the 
local order parameters are expected to follow the finite-size scaling form in the critical region around Tc (Eq. [3] ) 

[p(r;z,L;T)]9^«^ % ((T - T,)L^''') (39) 

For T < Tc, the convergence to finite- values [p{r; L — oo; T)]i in the L oo limit yields 

[p(r;L = oo;r)]9 = {Tc - Tf^"'' with /3(q) ^ iyf{q) (40) 

So the presence of a multifractal spectrum x{q) ^ qx{l) at criticality corresponds to non trivial exponents (3{q) ^ <7/3(l) 
for the powers of the local order parameters in the ordered phase. 



III. REMINDER ON WETTING AND POLAND-SCHERAGA TRANSITIONS 



A. Wetting and Poland-Scheraga models 

Wetting transitions are in some sense the simplest phase transitions, since they involve linear systems [1^. Let us 
consider a one-dimensional random walk (RW) of 2L steps, starting at z(0) = 0, with increments z(r+l) — z{r) = ±1. 
The random walk is constrained to remain in the upper half plane z > 0, but gains an adsorption energy if z{r) = 0. 
More precisely, the model is defined by the partition function 

^wettinq 

{2L) = ^ exp (3 ^ <ir5z2^fi (41) 

RW \ l<r<L J 

with inverse temperature (3 ~ l/T. In the pure case — there exists a continuous phase transition between a 
localized phase at low temperature, characterized by an extensive number of contacts at z = 0, and a delocalized 
phase at high temperature. 
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The Poland-Scheraga (PS) model of DNA denaturation [3g| is closely related to the wetting model. It describes 
the configuration of the two complementary chains as a sequence of bound segments and open loops. Each loop of 
length / has a polymeric entropic weight Af{l) ^ mV^'^j whereas each contact at position r has a Boltzmann weight 
g-/3er^ We assume that the two chains are bound at r = 1 and r = L. The partial partition hmction Zps{r) with 
bound ends then satisfies the simple recursion relation 

r-l 

Zps{r) - e-P^^' J2 -^(^ - ^')Zps{r') (42) 

r' = l 

The wetting model (|4T|) corresponds to a Poland-Scheraga model with loop exponent c = 3/2 (this exponent comes 
from the first return distribution of a one-dimensional random walk) . For DNA denaturation, the appropriate value 
of the loop exponent c has been the source of some debate. Gaussian loops in d = 3 dimensions are characterized by 
c = d/2 — 3/2. The role of self avoidance within a loop was taken into account by Fisher [33], and yields the bigger 
value c = dvsAW ^ 1.76, where vsAw is the SAW radius of gyration exp onent in rf = 3. More recently, Monte Carlo 
simulations of self avoiding walks jsl, [3§| and theoretical arguments [iol pointed towards a value c > 2. 

B. Disorder relevance as a function of the loop exponent c 

The Harris criterion concerning the stability of pure second order transitions with respect to disorder relies on the 
sign of the specific heat exponent 

2c - 3 

ctp = 2- vp = 43) 

c — 1 

Disorder is thus irrelevant for 1 < c < |, marginal for c = 3/2 relevant for | < c < 2. Poland-Scheraga models 
are thus particularly interesting to study disorder effects on pure phase transitions, since the parameter c allows to 
study, within a single model, the various cases of second order transition with respectively marginal/relevant disorder 
according to the Harris criterion, or first-order transition. From this point of view, it is reminiscent of the 2D Potts 
model, where the pure critical properties vary with the parameter q : the transition is second order for g < 4, the 
Ising case q — 2 corresponding to the marginal case of the Harris criterion, whereas the transition becomes first order 
for q > 4. The marginal case c = | has been studied for a long time [H, IH, [H, [H, [H, [H, |4^ and is of special 
interest since it corresponds to two-dimensional wetting as explained above. 

C. Numerical details 

In the following, we will study the multifractal properties of the local contact density 

p{r;i,L)=<S^^^o>,.L (44) 

representing the probability that the monomer r of the sample (i) of length L is on the interface z = at criticality 
T = Tc- We have chosen the same disorder distribution and parameters as in our previous work [49| . and we have used 
the same Fixman-Freire scheme to speed up calculations, as explained in details in [H, |4^. The results presented 
below have been obtained for the following sizes L and the corresponding number ns{L) of disordered samples 

= 16,32,64,128,256,512 (45) 



103 

ns{L) 



500,250,120,60,30,15 (46) 



IV. MULTIFRACTAL ANALYSIS OF THE WETTING TRANSITION WITH LOOP EXPONENT c = 1.75 



In this Section, we describe our results for the wetting transition with loop exponent c — 1.75 that corresponds to 
relevant disorder as explained above (Eq. 115)) . 
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FIG. 1: (Color online) Wetting with loop exponent c = 1.75 (a) Exponents x{q) governing the decay of the disorder averaged 
q-th powers of the local order parameter (Eq. I47|l (b) Generalized dimensions D{q) (Q) and D{q) (□) associated respectively 
to the typical values (Eq. [8]) and to disorder- averaged values (Eq. [T9)l of the generalized moments (Eq. [51) . 



A. Exponents x{q) and generalized dimensions D{q) and D{q) 

We shovif on Fig. [T] (a) the scaling dimensions x{q) governing the disorder-averaged moments of the local contact 
density (Eq. [3]) 

1 



< '^^..0 >'^^ (47) 

It is strongly non linear x{q) ^ <ix(X)- Some particular values are 

a;(q = 1) ~ 0.63 (48) 

x{q = 2) ~ 0.71 (49) 

In the large q limit, it saturates towards 

x{q H-oo) 1 (50) 



This point will be discussed in details in Section rVTl 

On Fig. [T](b), we show the generalized dimensions D{iq) and D{(i) associated respectively to the typical values (Eq. 
[5]) and to disorder-averaged values (Eq. [TO)) of the generalized moments (Eq. ^ 

Y^(r^L) = Y yi (51) 

In particular, the information dimension of Eq. 1101 is 

15(5 = 1) = i)(l) ~ 0.59 (52) 

and the correlation dimension of Eq. [Tl] is 

B{q = 2) - b{2) ~ 0.44 (53) 
For q large enough, the two exponents do not coincide anymore 13(g) 7^ D(q), as expected from the discussion on Eqs 

EH El El 
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FIG. 2: (Color online) Wetting with loop exponent c — 1.75 (a) Typical singularity spectrum f{a) (Eq. I14p : the maximum 
occurs at uq ~ 1.53 which represents the typical value. The minimal value is around Qmm — 0.36. (b) Corresponding curve 
a{q). 

B. Typical singularity spectrum /(a) 

To measure the typical singularity spectrum introduced in Eq. 1141 we have used the standard method based on 
g-measures of Ref. 51]. We show on Fig. [21(a) the curve f{a). The maximum corresponds to the typical exponent 
ao (Eq. [HI) 



atyp = CtQ = 1.53 

The curve is tangent to the diagonal a — f{a) at the point (Eq. [T5)l 

ai = /(ai) = i?(l) ~ 0.59 

The minimal value corresponds to 

amin = D{q = +oo) ~ 0.36 

and the maximal value 

amax = D{q = -oo) ~ 2.22 



(54) 
(55) 
(56) 
(57) 



On Fig. [2](b), we show the corresponding curve a{q) representing the dominant exponent a that contribute to the 
g-generalized moment (Eq. I15p . 



C. Disorder-averaged singularity spectrum f{a) 

In contrast to the method of 51] that allows to measure numerically the typical spectrum, we are not aware of an 
efficient method to measure the disorder-averaged singularity spectrum f{a). We have thus measured the probability 
distributions Hl {a) of the rescaled weights 



In w(r; i, L) 
hTE 



(58) 



in analogy with similar numerical measures of the multifractal spectrum from the statistics of correlation function 
in disordered Potts models [l^ [13] • Our results presented on Fig. [3] show that the convergence towards the typical 
spectrum /(a) in the positive region a > is extremely slow. In particular, the convergence of the most probable 
exponent amp{L) towards the typical value ao is extremely slow, of order l/ln(L). 
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In w 
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In L 

FIG. 3: (Color online) Wetting with loop exponent c = 1.75 : the histograms Hl{ol) of the rescaled variable a — _ lfli£lljj2ii) 
shown for 2.10"^ < ^ < 512.10^ converge extremely slowly (as l/ln(L)) towards the typical singularity spectrum /(a) (O) 
measured with the method of Ref . [sJl • 




FIG. 4: (Color online) Wetting with loop exponent c = 1.75 (a) histogram over the samples of the information entropy s{i,L) 
defined in Eq. [TO]for sizes 10^ < L < 512.10^ : the averaged value grows as s{i,L) ~ -D(l) InL (b) histogram over the samples 
of lnF2(i,i) (Eq. |6]) for sizes 10^ < L < 512.10^ : the averaged value behaves as \nY2{i,L) ~ -Z)(2)lnL. 

D. Histograms of the 'information' entropy and of I2 

We show on Fig. [4] (a) the histogram over the samples (i) of the 'information' entropy s{i, L) defined in Eq. [10] : 
as L grows, the averaged value grows logarithmically (Eq. llOp 

s{i,L) ~ £1(1) InL (59) 

whereas the width converges towards a constant value. 

Similarly, we show on Fig. |4] (b) the histogram over the samples (i) of In Y2 (?, L) (Eq. [6] ) : as L grows, the averaged 
value grows logarithmically (Eq. [TT|) 



\nY2{i,L) ~ -D{2)lnL 
whereas the width converges towards a constant value. 



(60) 
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FIG. 5: (Color online) Wetting with loop exponent c = 1.75 : Wiseman-Domany lack of self-averaging at criticality for the 
spatial average p{i, L) of the order parameter (Eq. [l| (a) histogram over the samples of the variable u = p{i, L) / p{i, L) for sizes 
16.10^ < L< 512.10^ : u remains a random variable of order 0(1) in the limit L oo. (b) same data in logarithmic scale. 



E. Statistics of the spatial averaged order parameter over the samples 



To study the Wiseman-Domany lack of self-averaging at criticality for the spatial average L) of order parameter 
(Eq. [28]) , we have computed the probability distribution Gl [u) of the ratio 

u=P=^ (61) 
p{i,L) 



The results for various L presented on Fig[S]show that u remains a random variable of order 0(1) in the limit L — )■ oo. 



V. MULTIFRACTAL ANALYSIS OF THE WETTING TRANSITION WITH LOOP EXPONENT c = 1.5 



In this Section, we describe our results for the wetting transition with loop exponent c = 1.5 that corresponds to 
marginal disorder as explained above (Eq. |43)) . 

A. Exponents x{q) and generalized dimensions D{q) and D{q) 

We show on Fig. [S] (a) the scaling dimensions x{q) governing the disorder-averaged moments of the local contact 
density (Eq. [3]) 

< <5..,o >' « (62) 
It is strongly non linear x{q) ^ ^.^(l). Some particular values are 

x{q = 1) ~ 0.52 (63) 
x{q = 2) ~ 0.79 (64) 

In the large q limit, it saturates towards 

x{q +oo) — > 1 (65) 

(see the discussion of Section rvT|) . 
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FIG. 6: (Color online) Wetting with loop exponent c — 1.5 (a) Exponents x{q) governing the decay of the disorder averaged 
g-th powers of the local order parameter (Eq (I62|l (b) Generalized dimensions D{q) (Q) and D{q) (□) associated respectively 
to the typical values (Eq. [8]) and to disorder- averaged values (Eq. [T9)l of the generalized moments (Eq. |66)) . 



On Fig. [6](b) we show the generalized dimensions D{q) and D{q) associated respectively to the typical values (Eq. 
[8]) and to disorder-averaged values (Eq. [T9|) of the generalized moments (Eq. [6|) 

YS.L) = (66) 

(5;<^.,»>) 

In particular, the 'information dimension of Eq. [10] is 

£i(g = 1) = i)(l) ~ 0.87 (67) 

and the correlation dimension of Eq. [11] is 

D{q = 2) - b{2) ~ 0.74 (68) 



B. Typical singularity spectrum /(a) 

We show on Fig. [7] (a) the curve f{a) obtained via the standard method of Ref. jSli] The maximum corresponds 
to the typical exponent ao (Eq. [T7|) 

atyp = 0.0 = 1-125 (69) 
The curve is tangent to the diagonal a — f[a.) at the point (Eq. [18]) 

ai =/(ai) =i?(l) - 0.87 (70) 

The minimal value corresponds to 

amin = D{q = +oo) ~ 0.48 (71) 

and the maximal value 

amaa; -OO) ~ 1.53 (72) 

On Fig. [3(b), we show the corresponding curve a{q) representing the dominant exponent a that contribute to the 
g-generalized moment (Eq. [T5]l . 
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FIG. 7: (Color online) Wetting with loop exponent c = 1.5 (a) Typical singularity spectrum f{a) (Eq. I14p : the maximum 
occurs at ao ~ 1.125 which represents the typical value. The minimal value is around — 0.48. (b) Corresponding curve 

a{q). 
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FIG. 8; (Color online) Wetting with loop exponent c — 1.5 : the histograms Hl(oi) of the rescaled variable a — _]}12!^^j1i£l 
shown for 2.10'^ < < 512.10^ converge extremely slowly (as l/ln(L)) towards the typical singularity spectrum /(a) (O) 
measured with the method of Ref. 51]. 



C. Disorder-averaged singularity spectrum /(a) 
We show on Fig. [8] the histograms i?L(a) of 

lnw(r;z,£) 
" = 

for various sizes L and compare with the the typical spectrum / (a) obtained via the method of Ref. '5l'|. As previously 
mentioned for the corresponding Figure [3] concerning the case c — 1.75, the convergence is extremely slow, of order 
l/h^(L). 



D. Histograms of the 'information' entropy and of Y2 

We show on Fig. [9] (a) the histogram over the samples (i) of the 'information' entropy s(i, V) defined in Eq. [10] 
as L grows, the averaged value grows logarithmically (Eq. [TOl) 



s{i,L) ~ £1(1) InL 



(74) 



14 



L=1000 



3 



I A A A A 



10 



(a) 



U512000 



12 

s(i,L) 




FIG. 9: (Color online) Wetting with loop exponent c = 1.5 (a) histogram over the samples of the information entropy s{i,L) 
defined in Eq. [TO]for sizes 10^ < L < 512.10^ : the averaged value grows as s{i,L) ~ D{1) InL (b) histogram over the samples 
of lny2(i, L) (Eq. |6]) for sizes 10^ < L < 512.10^ : the averaged value behaves as lnY2ii,L) ~ -D{2) InL. 



whereas the width again converges towards a constant value. 

Similarly, we show on Fig. [9] (b) the histogram over the samples (i) of In Y2 {i, L) (Eq. [6] ) : as L grows, the averaged 
value grows logarithmically (Eq. [Tl|) 



liiY2ii,L) ~ -D{2)liiL 
whereas the width converges towards a constant value. 



(75) 



E. Statistics of the spatial averaged order parameter over the samples 




FIG. 10: (Color online) Wetting with loop exponent c = 1.5 : Wiseman-Domany lack of self-averaging at criticality for the 
spatial average p{i, L) of the order parameter (Eq. [l| (a) histogram over the samples of the variable u = p{i, L)/p{i, L) for sizes 
16. 10'^ < L < 512.10^ : u remains a random variable of order 0(1) in the limit L 00. (b) same data in logarithmic scale. 



To study the Wiseman-Domany lack of self-averaging at criticality for the spatial average p(i, L) of order parameter 
(Eq. [55]), we have computed the probability distribution Gl{u) of the ratio 

u = ^ (76) 
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The results for various L are shown on Fig 1101 show that u remains a random variable of order 0(1) in the limit 
L —t oo. 

VI. INFLUENCE OF BOUNDARY CONDITIONS 

A. Value x[q oo) and minimal value Omin of the multifractal spectrum 

In our numerical studies for c — 1.75 and c = 1.5 presented above, we have found in both cases that the exponents 
x{q) governing the decay of the powers of the local order parameter (Eq[3]) saturate to 1 for large q (EgslSOl and [65|) 



x{q oo) = 1 



(77) 



This indicate that moments of high order g ^ oo are dominated by a finite number of points (of density of order l/L) 
having a finite order parameter 0(1). So the minimal ymin actually saturates the bound of Eq. 1381 



Urain — ^min d x(^\^ — 

As a consequence, the minimal value amin of the multifractal spectrum is simply given by 

O^min d ^(1) 



(78) 



(79) 



This simple relation is satisfied with the above numerical results both for c = 1.75 (with x{q = 1) ~ 0.63 and 
amin = D(q = +oo) ~ 0.36) and for c = 1.5 ( with x{q = 1) ~ 0.52 and amin = D{q ~ +oo) ~ 0.48 ) 

To determine whether these properties are linked to the bound-bound boundary conditions used, or are more 
general, we have studied other boundary conditions as we now explain. 



B. Comparison between fixed and free boundary conditions 



f (a) 





FIG. 11: (Color online) Typical singularity spectra f{a) for bound-bound (Q) or free-free (□) boundary conditions (a) wetting 
with loop exponent c = 1.75 (b) wetting with loop exponent c — 1.5. 



In the previous Sections IIVI and |Vl we have presented numerical results concerning the bound-bound boundary 
conditions, where the polymer is attached to the interface z = at both ends 



z(r = 1) = z{r ^ L) 



(80) 



Since the wetting models cannot be directly defined with completely free boundary conditions (since the space above 
the interface is infinite), we have still considered the bound-bound boundary conditions of Eq. [50] but we have 
measured the multifractal spectrum using only the 'bulk monomers' satisfying 



L 3L 
— <r < — 
4 - - 4 



(81) 
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This procedures aims to simulate 'free' boundary conditions at r = j and at r = ^ for a sample of size ^. The 
singularity spectra measured with these two types of boundary conditions are shown on Fig. [TT] (a) and (b) for 
the cases c = 1.75 and c = 1.5 respectively. These singularity spectra coincide within our numerical accuracy (the 
difference on the right half of Fig. [TT] (a) corresponds to the negative moments q < Q dominated by the smallest 
weights whose statistics is more difficult to measure precisely). 



C. Griffiths ordered clusters 



In particular, the fact that the minimal value ttmin remains the same for fixed and 'free' boundary conditions shows 
that the finite number of points having a finite order parameter are not confined to the boundaries but also exist in 
the bulk at criticality. This is related to the finite probability of finite ordered clusters in diluted disordered systems 
below the pure critical temperature, which have been much studied in the context of Griffiths singularities [52| . 

Our conclusion for the wetting transition and more generally for disordered systems at criticality, is that the finite 
probability of Griffiths ordered clusters determines the asymptotic value 

xiq oo) ^ d (82) 

of the exponents x{q) governing the decay of the powers of the local order parameter. Accordingly, the minimal Umin 
saturates the bound of Eq. [38l 

Umin = amm - d + x{l) = (83) 

and the minimal value amin of the multifractal spectrum is simply related to the exponent a;(l) governing the decay 
of the order parameter 

armn^d~x{l) (84) 



VII. SUMMARY AND CONCLUSIONS 



In this paper, we have studied in detail the multifractal statistics of the local order parameter at random wetting 
critical points, for two values of the loop exponents c, namely c = 1.75 and c — 1.5. For these models where large 
sizes up to L = 512.10'^ can be probed with a good statist ics ov er the samples, we have numerically measured 

(i) the 'Ludwig exponents' x{q) that govern the moments p'^{r) of the local order parameter p(r) 

(ii) the generalized dimensions D(q) and D{q) associated to typical and disorder-averaged values of the moments 
Yq = w''{r) of the multifractal measure w{r) = p{r)/{J2r' 

(iii) the corresponding typical singularity spectrum /(a). 

We have also discussed the relations between this multifractal statistics and the Wiseman-Domany lack of self- 
averaging at criticality. Finally, we have argued that the presence of finite Griffiths ordered clusters at Tc determines 
the asymptotic value of the Ludwig exponent x{q oo) = d and the minimal value amin — D{q oo) = d — x{l) 
of the multifractal spectrum. We have checked that these relations are well satisfied in our numerical results for the 
wetting transitions (where d = 1), not only for bound-bound boundary conditions but also for 'free-free' boundary 
conditions. 

However finite Griffiths ordered clusters occur above critical points in disordered systems independently of the 
relevance/irrelevance of the disorder from the Harris criterion ^52.], which concerns coarse-grained properties at weak 
disorder. Our conclusion is thus that the multifractal statistics of the local order parameter should be non-trivial for 
any critical point with frozen disorder, since it probes the heterogeneities at all scales. This would explain why our 
results concerning the marginal disorder case c = 1.5 are qualitatively similar to our results for the relevant disorder 
case c = 1.75. 
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